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Abstract 

For a slowly rotating non magnetized white dwarf the accretion disk extends all the way to the star. Here the matter impacts 
and spreads towards the poles as new matter continuously piles up behind it. We have solved the 3d compressible Navier-Stokes 
equations on an axisymmetric grid to determine the structure of this boundary layer for different viscosities corresponding to 
different accretion rates. The high viscosity cases show a spreading BL which sets off a gravity wave in the surface matter. 
The accretion flow moves supersonically over the cusp making it susceptible to the rapid development of gravity wave and/or 
Kelvin-Helmholtz instabilities. This BL is optically thick and extends more than 30 degrees to either side of the disk plane 
after 3/4 of a Keplerian rotation period (tx=19s). The low viscosity cases also show a spreading BL, but here the accretion 
flow does not set off gravity waves and it is optically thin. 
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1 Introduction 

Cataclysmic variables (CV) is an interesting class of 

■ close binary stars comprising a hot white dwarf (WD) 
, and a relatively lower mass red dwarf star filling its 

■ Roche lobe [36]. In such systems hydrogen-rich matter 
\ from the red dwarf exits through the inner Lagrange 
. point and flows towards the white dwarf, where the mat- 
' ter proceeds to form an accretion disk around the white 
, dwarf due to the excess angular momentum originating 

■ from the orbital motion of the binary [28] . 

■ In the disk, turbulence [33] and magnetic fields [2] dis- 
sipate potential energy and cause the matter to lose an- 
gular momentum. As a result the disk matter slowly 
falls inwards towards the white dwarf [29]. If the WD 
is non-magnetic, the accretion disk extends all the way 
to the surface. Only half of the entire accretion lumi- 
nosity, Lace = GM^,m/ R^,, is emitted in the disk, since 
the matter is still moving at a roughly Keplerian veloc- 
ity, vk ~ i/GM^/i?,, before it is ultimately accreted 
[22]. The kinetic energy must therefore be converted into 
other forms of energy in order for the matter to come 
into co-rotation with the surface of the WD [29] . This 
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occurs in a comparably small (on the size of a disk scale 
height) boundary layer (BL) near the star. 

Since the BL is much smaller than the disk itself and each 
radiates equal amounts of energy, the disk emits mainly 
in the optical and near UV, whereas the BL emits in the 
far UV and in soft X-rays. However, since the kinetic 
energy may also be converted into winds, WD rotation 
or heating, the details of the emitted spectrum depends 
on the detailed structure of the BL. Therefore the struc- 
ture of the BL is important to understand both to make 
the observational interpretations between the different 
components of the spectrum which is important in the 
study of dwarf novae and in the general interpretation 
of observational data, but also because it determines 
how the accreted matter ultimately distributes itself in 
the envelope of the WD, which is important for classical 
novae. Before such a detailed analysis of the radiative 
transfer in conjunction with the hydrodynamics can be 
undertaken, it is important to understand the dynam- 
ical processes the lead to the formation of the BL. We 
focus on such a dynamical study in this paper leaving a 
radiation-hydrodynamical study for later. 

Dwarf novae are a subclass of CVs in which a thermal 
instability in the accretion disk causes an increase in the 
matter transfer rate through the disk and thus an in- 
crease in the rate of gravitational energy release [9] . This 
in turn causes a 2-5 mag increase in the disk luminos- 
ity and consequently in the BL luminosity. The BL may 
react in various ways due to enhanced accretion and de- 
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tailed studies are required to determine its response [13] . 

Another CV subclass is the classical nova. A classical 

nova obtains, when enough hydrogen-rich rich material 
has been accumulated in the envelope to reach the igni- 
tion point at the e-degenerate base of the envelope. The 
resulting thermonuclear runaway ejects part or all of the 
envelope. The initial CNO composition of the burning 
material should be strongly enhanced compared to the 
accreted material to account for composition of the ob- 
served ejecta [34]. Several mechanisms have been sug- 
gested for this CNO enhancement (see [17] and refer- 
ences therein). They can be roughly be divided into pre- 
burst mixing between accreted material and the under- 
lying CO rich WD [18,23,32,1] by accretion driving in- 
stabilities or mixing with the underlying material when 
the convective zone of the thermonuclear runaway ex- 
tends deep enough to dredge up CO material [34,15]. 

The BL has typically been solved in a model where the 
averaging of the vertical structure which is employed 
in accretion disk studies [33] has been extrapolated 
to the surface. This reduces the description of the BL 
to a one dimensional problem in the radial direction 
[29,24,27,10]. However, this assumption is not neces- 
sarily correct since the BL may spread out [12]. These 
calculations were complimented with an analytic treat- 
ment of the meridional direction which showed that BL 
could spread and cover a significant part of the WD [26]. 

Recognizing the multi-dimensionality of the problem 
several authors have attacked the problem directly us- 
ing numerically methods culminating in simulations 
which included flux-limited radiative transport as well 
as different viscosities [30,21,19,20]. However, as these 
simulations predated those of [26] their setup did not 
include sufficient resolution in the radial and meridional 
planes to capture the dynamics and structure of the BL. 

In this paper we set out to calculate the dynamical struc- 
ture of the BL with sufficient numerical resolution to 
capture the dynamical evolution of the accretion flow 
and its interaction with the stellar surface (see [13]). Our 
model is presented in section 2 and results for a = 0.1, 
a = 0.03, a = 0.01, a = 0.005, and a = 0.001 are given 
in section 3. 

2 Computational model 

The source of the angular momentum transport in the 
BL, which is responsible for transporting matter through 
the disk, is a combination of magnetic fields and tur- 
bulence. However, an a priori prescription, in particu- 
lar in the BL, is a source of disagreement (see [27] and 
references therein) . Instead the efficiency of the angular 
momentum transport can be parametrized with a coef- 
ficent, a [33]. Describing the angular momentum trans- 
port with a simple shear coefBenct means that the dy- 



namics follows the Navier-Stokes equations. Here the 
Navier-Stokes equations as given by [25] are solved in 
spherical coordinates {r,6,<p) on an axisymmetric mesh 
with 384 ratioed zones in the radial direction and 128 
ratioed zones in the meridional range spanning to 30 
degrees from the disk plane same as [13]. Symmetry 
across the disk plane is assumed. This allows a doubling 
of the zone resolution at the same CPU-cost. 

This allows us to resolve a pressure scale height in the 
radial direction and a disk scale height in the meridional 
direction for a stellar temperature of T=300000K and 
a disk temperature of T=1000000K. We assume a thin 
halo which consist of the same material as the disk (and 
the star). The halo provides a pressure that balances the 
pressure at the upper boundary of the disk, keeping the 
upper boundary of the disk in equilibrium. The disk and 
halo model are described in [5] . 

The spatially and temporally second order agorithms 
in RIEMAMM have been described in [31,3,4,7,8,5] and 

use many ideas from higher order WENO schemes (see 
[16,6]) to reduce dissipation. The matter in the model is 
subject to the central gravitational field of the underly- 
ing WD, which has a radius and mass of i?* = 9x10* cm 
and M* = 0.6Mq respectively - same as [26]. The model 
uses an ideal gas (7 = 5/3) of a fully ionized solar com- 
position (/X = 0.62g/mole) and assumes no radiation 
transport. 

3 Results &: Discussion 

In the next few subsections, we focus on several impor- 
tant aspects of the boundary dynamics. In section 3.1, 
we focus on the evolution of angular momentum in the 
boundary layer. In section 3.2, we discuss accretion 
based instabilities which may obtain and in section 3.3 
we relate the computations to observations. 

3. 1 Angular momentum evolution 

As the simulations start, the shear force transfers angu- 
lar momentum between differentially rotating disk an- 
nuli so that matter can move inwards and accrete on the 
star at the footpoint of the disk. Angular momentum 
transfer between the disk material and the stellar surface 
is also necessary so matter can move to higher latitudes. 
Otherwise the centrifugal barrier prevents matter from 
leaving the footpoint of the disk. It is this transfer which 
spins up the existing surface layers of the star. Angular 
momentum is also directly advected onto the star, be- 
cause the orbiting disk material eventually accretes and 
forms the new surface. 

Fig. 1 shows the resulting specific angular momentum 
after 3/4 of a Keplerian evolution for the a = 0.1 case, 
where fig. 2 shows it for the a = 0.001 case. The halo is 
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spun up as is the underlying star at the footpoint, where 
the disk connects. At higher lattitudes the spreading disk 
material is less dense and thus takes longer amounts of 
time to drag the star's surface around with it. 
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Fig. 1. A contour plot of the specific angular momentum, 
Q = v^/r, for the a — 0.1 case after 3/4 of a Keplerian 
rotation period. 
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Fig. 2. A contour plot of the specific angular momentum for 
the Q — 0.001 case after 3/4 of a Keplerian rotation period. 

Since degenerate matter is too "slippery", turbulence 
and/or threaded magnetic fields might be the two most 
prominent ways to transfer a significant amount of an- 
gular momentum to the deeper layers of the WD [11], 
where the relative motion drives the mixing between the 



accreted surface material and the deeper layers. A su- 
personic component in the toroidal direction over most 
of the BL obtains for all values of a (see fig. 3). Yet 
only the a = 0.1 case has a supersonic component in 
the poloidal directions, whereas the poloidal flows for 
a = 0.03, a = 0.01, a = 0.005, and a = 0.001 remain 
subsonic (see figs. 1 and 2) [13] . The supersonic velocities 
mean that the flows are susceptible to the rapid develop- 
ment of gravity wave and/or Kelvin Helmholtz instabili- 
ties. Such instabilities would develop self-consistently, if 
the interface had been given substantially more zones in 
which to form these instabilities. Due to our limited al- 
beit high zone resolution, wc can only provide the veloci- 
ties which provide the driving mechanism for such insta- 
bilities. The calculation of the turbulent mixing which 
obtains from these instabilities must therefore be calcu- 
lated using other models [18,23,32,1]. 



Fig. 3 shows the torodial velocity at the meridional cen- 
terline as a function of radius for the simulated alpha- 
viscoties. The cases of a = 0.01, a = 0.005, and a = 
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Fig. 3. Toroidal (<j!>) velocity at the meridional centerline of 
the disk shows the extent of the boundary layer as well as 
the sound speed after 3/4 of a Keplerian rotation period. 
Note that asymptotes to a higher value than Cs,a. 

0.001 do not show any significant spin-up of the star 
after 3/4 of a Keplerian rotation, whence there is not 
significant shear connection between the disk and the 
stellar surface. Therefore disk matter retains its Keple- 
rian motion much closer to the star which means that 
even though the accretion rates are smaller, significant 
amounts of the dissipation can still be generated in the 
BL close to the star. 

3.2 Formation of gravity waves 

For the a = 0.1 case, the edge of the boundary layer 
{dn/dr = 0, see [14]) is located at 1.06i?*. At i? = 
1.03i?*, is about 92% of u^. This means that only 84% 
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of the gravitational pull is supported by the centrifugal 
force while the rest is supported by pressure. The pres- 




Fig. 4. A contour plot of the pressure for the a — 0.1 case 
after 3/4 of a Keplerian rotation period shown with a loga- 
rithmic scale. 

sure for a = 0.1 is illustrated in fig. 4 and comes from 
the build up of density due to the high accretion rate 
facilitated by the high viscosity which drags down mate- 
rial from the disk and also keeps it from moving rapidly 
towards the poles once it makes contact with the WD 
surface. This results in a dense band at the foot point of 
the disk which causes a gravity wave of surface matter to 
spread towards the poles. The matter inflow is, however, 
sufficiently large to cause the disk material to overflow 
the gravity wave supersonically as described above and 
shown in fig. 1. Over time the thick band will spin up 
the surface and thus prevent shear forces from accreting 
matter since the differential rotation decreases as the BL 
is spun up. 

The propensity for the creation of gravity waves and 
their magnitude decreases rapidly with lower values of 
a. Even for a = 0.03, the effect is only marginal and for 
lower values it is no longer present. The effect is thus 
only present during dwarf novae in the outburst stage. 
Furthermore the gravity wave might be transient as the 
surface adjusts to a fluctuations in the accretion rate. 
For lower values of a, the matter accretes slowly inwards 
and spreads towards the poles in a uniformly thick layer. 

3.3 Dissipative heating and optical depth of the BL 

Observationally, the high density of the spreading BL 
leads to optically thick conditions for a = 0.1 and a = 
0.03 in a band extending at least 30 degrees away from 
the disk's centerline. The band should therefore provide 



a continuum component to the observed spectrum. Here 
we use the criterion that — Kpdiskdr > 2/3, where 
K = 0.34 gcm"-^ is the Thomson scattering opacity of 
a fully ionized solar composition. This constitutes the 
minimum amount of scattering and thus provides a lower 
bound of the opaqueness of the matter. Since the density 
is lower for lower values of a, the BL is optically thin 
more than 15 degrees away from the disk itself for a = 
0.01, a = 0.005, and a = 0.001. 

The local rate of energy dissipation depends on the di- 
vergence of the angular velocity [25] , so flg. 3 also indi- 
rectly shows where heat is released and thus where ra- 
diation is emitted. Fig. 5 shows that the supersonic im- 
pact of the accretion flow following the cusp dissipates 
the most energy. However, energy is also dissipated at 
the interface between the BL and the surface. This is 
better seen in fig. 6 which shows that the toroidially ro- 
tating matter dissipates rotational kinetic energy along 
its entire interface with the surface as it slides against it. 
This may explain the optically thin BL which may ex- 
plain the line emission as seen in observations of U Gem 
during the quiescent phase [13]. 




Fig. 5. A logarithmic contour plot of the dissipation rate for 
the a — 0.1 case after 3/4 of a Keplerian rotation period. 



4 Conclusion 

We have numerically simulated the structural dynamics 
of the BL for an accreting white dwarf surrounded by 
an a-disk for different values of a. The structure and 
dynamics of the BL is important, because it determines 
the specifics of the radiated spectrum emitted from this 
region which may account for up to half of the total 
energy released in the accretion process. For high values 
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Fig. 6. A logarithmic contour plot of the dissipation rate for 
the Q — 0.001 case after 3/4 of a Keplerian rotation period. 
Notice that the color scale is different for fig. 5. 

of a, the BL is optically thick and extends more than 
30 degrees to either side of the disk plane after 3/4 of a 
Keplerian rotation period (i/f=19s). 

The dynamics of the accretion process in the bound- 
ary layer is also important, because it determines the 
mixing ratio between the accreted matter and the un- 
derlying surface matter and thus determines the initial 
composition which is important for nucleosynthesis in 
classical novae. Here the simulations show that high val- 
ues of a result in a spreading BL which sets off gravity 
waves in the surface matter, where the accretion flow 
moves supersonically over the cusp making it suscep- 
tible to the rapid development of gravity wave and/or 
Kelvin-Helmholtz instabilities. The low viscosity cases 
also show a spreading BL, but here the accretion flow 
does not set off gravity waves and it is optically thin. 

The cooling time for the plasma is very short for all the 
simulations [35] , so in an upcoming paper by Balsara & 
Fisker we repeat the calculations without a viscous dis- 
sipation in the energy equation. This corresponds to an 
instant steady state cooling of the layer. Incorporating a 
flux limited diffusion radiation transport will be defered 
to a later paper. 
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